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Context. PSRB1259-63/LS2883 is a binary system where a 48 ms pulsar orbits a massive Be star with a highly eccentric orbit 
(e = 0.87) with a period of 3.4 years. The system exhibits variable, non-thermal radiation visible from radio to very high energies 
(VHE) around periastron passage. This radiation is thought to originate from particles accelerated in the shock region between the 
pulsar wind (PW) and stellar outflows. 

Aims. The consistency of the H.E.S.S. data with the inverse Compton (IC) scenario is studied in the context of dominant orbital phase 
dependent adiabatic losses. 

Methods. The dependence of the observed TeV flux with the separation distance is analyzed. Model calculations based on IC scattering 
of shock accelerated PW electrons and UV photons are performed. Different non-radiative cooling profiles are suggested for the 
primary particle population to account for the variable TeV flux. 

Results. The TeV fluxes obtained with H.E.S.S. in the years 2004 and 2007 seem to be only dependent on the binary separation. The 
presented results hint at a peculiar non-radiative cooling profile around periastron dominating the VHE emission in PSRB1259-63. 
The location of the stellar disc derived from this non-radiative cooling profile is in good agreement with that inferred from radio 
observations. 
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1. Introduction 

The binary system PSRB1259-63/LS2883 comprises a pulsar 
with a period of 48 ms and a massive star of spectral type 
B2Ve featuring a win d with a polar and equatorial component 
dJohnston et al.lll992l) . The latter slower and denser wind consti- 
tutes a circumstellar disc. The interaction of the pulsar with the 
dense environment (radiation and/or mat ter) in the binary system 
may give rise to non-thermal processes. iTavani & Aronsl (1 1997b 
have discussed various plausible scenarios for the production of 
non-thermal emission in the system. The following will focus on 
the no n-thermal processes that o ccur in the shocked pulsar wind 
(PW) dKennel & Coroniti|[l984T) . Indeed, the supersonic PW is 
terminated at a strong shock wave, and at this shock particles 
may be accelerated up to very high energies (VHE). These ac- 
celerated particles could, in case of a leptonic PW, upscatter pho- 
tons from the companion star to VH E via inverse Com pton (IC) 
processes producing VHE emission (Tavani et afll994l) . 

In the framework of this scenario, iKirk et all d 1999b have 
predicted the level and orbital phase dependence of the y-ray 
flux. Two regimes of non-thermal particle cooling were con- 
sidered: (i) negligible non-radiative (adiabatic) energy losses 
(i.e. dominance of radiative losses); and (ii) dominant (or- 
bital phase independent) adiabatic losses. However, the obser- 
yations of PSR B 1259 -63 with the H.E.S.S. telescope array 
dAharonian et aj] l2005h found a lightcurve that significa ntly de - 
viates from the two lightcurves obtained by IKirk et al.l (Q999). 
This discrepancy supported hadronic scenarios for the produc- 
tion of VHE emission in the system. In the framework of this 
scenario, VHE protons originating from the pulsar interact with 
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the dense outflow from the companion star, producing VHE 
emission. A characteristic feature of this scenario are two sharp 
maxima in the VHE emission related to the pulsar passage 
throu gh the dense equatorial outflow of the star (Kawachi et al. 
2004). The importance of the dense circumstellar disc seems also 
to be indicated by the lac k of pulsed radio emi ssion shortly be- 
fore and after periastron (Johnston et al1 l2005l) . Moreover, cor- 
relations between the X-ray and TeV fluxes seemed to support 
the ha dronic scenario for the VHE radiation dChernvakova et ail 
2006), since in the case of proton-proton interactions not only 
VHE y-rays are produced, but also secondary leptons. These 
secondary particles may effectively emit through the IC and syn- 
chrotron channels providing the X-ra y and radio emission cor - 
related with the VHE energy band dChernvakova et aD 12006). 
Thus, the TeV flux from PSRB1259-63 should be related to 
the rate of primary interactions, i.e. to the location of the cir- 
cumstellar disc on the pulsar orbit. However, it should be noted 
that the reconstruction of the disc properties from the eclipse 
of radio pulses leads to somewhat different values for the geo- 
metrical parameters of the stellar disc than those required in the 
hadronic scenario (Bogomazov 2005). Moreover, H.E.S.S. data 
taken around the 2007 periastron covering orbital phases differ- 
ent to the ones in 2004 seem to contradict such a hadronic disc 
scenario since the TeV flux does not fo llow a simple disc den - 
sity profile as indicated by 2004 data dAharonian et al.l 12009). 
Thus, the currently obtained lightcurves apparently do neither 
support leptonic nor hadronic scenarios for the production of 
VHE photons in PSRB1259-63. This discrepancy suggests the 
existence of additional effects which have to be considered in 
the modelling. In the case of the hadronic scenario, the sim- 
plest assumption is the orbit-to-orbit variability of the disc lo- 
cation. This hypothesis may be supported by the orbit-to-orbit 
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variability of the p ulsed radio emission detected from the system 
dWang et al.ll2004l) but seems unlikely since there are no signif- 
icant inter-orbital variations of X-ray and unpu lsed radio fluxes 
dChernvakova et al.ll2009t I Johnston et aDl2005h . Clarifying this 
contradictory behavior with respect to a hadronic disc scenario 
requires a dedicated study. In what follows, a possible origin of 
the complicated behavior of VHE fluxes in the framework of a 
leptonic scenario will be discussed. 

The IC TeV energy flux is expected either to peak at peri- 
astron, where the target photon field density is maximal for 
the standoff point, i.e. at the location of the pulsar-star wind 
shock, or to show a rather smo oth orbital pha se dependence 
in case of a saturation regime dKirket al.lll999h . This is, how- 
ever, in opposition to what was measured by H.E.S.S. The 
measured VHE lightcurve clearly shows a drop in photon flux 
towards periastron. Moreover, the orbital phas e dependence of 
the X-ray lightcurve (Chernyakova et al. 2006) does not fit th e 
shape expected in a simple leptonic model dKirk et alj [l999). 
All these discrepancies a sk for a more sophisticated modelling. 
Khangulv an et al.l (|2007) modelled the VHE data obtained with 
H.E.S.S. in 2004 with synchrotron/IC scenarios that were ac- 
counting for (i) a change of the acceleration rate due to en- 
hanced losses close to the periastron passage, (ii) dominant or- 
bital phase dependent adiabatic losses, (iii) the orbital phase de- 
pendent change of the pulsar wind bulk Lorentz factor due to the 
intera ction with stellar photons. More recently, Uchivama et alj 
(2009) modelled the X-ray data obtained with the Suzaku X-ray 
satellite with a synchrotron and IC radiation mod el that included 
orbita l phase dependent adiabatic losses. Finally, Takata & Taam 
(2009) suggested that the observed variability of non-thermal 
fluxes is related to the orbital phase dependence of (i) the pul- 
sar wind bulk Lorentz factor and magnetization parameter, (ii) 
the magnetization parameter and the spectral power law index of 
shock-accelerated particles. In this last approach, non-radiative 
losses were not accounted for. 

The introduction of orbital phase dependent non-radiative losses 
accounts for the following physical effects most likely being 
present in PSRB 1259-63: (i) the post-shock flow propagating 
in the co nfined region formed by the PW and the shocked stel- 
lar wind (Bogovalo v et alJ 2008) should suffer significant adia- 
batic losses dKhangulvan et alj |2008): (ii) the loss rate depends 
on the hydrodynamics of the interaction and may thus natu- 
rally depend on the orbital separation distance and the density of 
the stellar win d. Moreover, since the pulsar wind is expected to 
be anisotropic (Bogovalov & Khangoulvan 2002), the adiabatic 
loss rate could have a rather complicated or bital phase depen- 
dence . Hence, in the following the approach of Khangulvan et al. 
(I2007h is used, taking into consideration the data obtained with 
H.E.S.S. in 2007. 

This paper is organized as follows. First of all, a phenomenologi- 
cal discussion of the PSR B 1259-63 lightcurve both in the y-ray 
and X-ray bands is presented in Sec. [2] A comparison of the ex- 
pected IC TeV energy flux with the H.E.S.S. data allows an infer- 
ence of the non-radiative cooling profile. This ansatz yields re- 
sults for the prediction of TeV photons from PSRB 1259-63 that 
describe the peculiar dro p at periastron with re asonable accuracy 
dKhangulvan et al.ll2007l) . In fheir[2007 paper. iKhangulvan et al.l 
(2007j) used the H.E.S.S. 2004 data in order to calculate the non- 
radiative cooling coefficients for the PSRB 1259-63 orbit. In 
Sec. [3] of this article their analysis is redone with the addition 
of VHE data from the 2007 periastron passage. Moreover, three 
different possible cooling profiles are compared. 



2. The lightcurve of PSR B1 259-63 as seen in 
X-rays and y-rays 



While y-ray binaries, such as LS5039 and LSI+61303, have 
been observed over several orbits and thus could be confirmed 
as pe riodical VHE emitters (Aha ronian et alj|2006t lAlbert et al.l 
2009) this was not possible for PSRB1259-63 so far because 
of the rather long orbital period of 3.4 years. The overlap be- 
tween the H.E.S.S. datasets of 2004 and 2007 is marginal and 
thus inconclusive in this regard (see top panel in Fig. [TJ. It is 
thus up to now unknown whether the TeV lightcurve of this ob- 
ject shows periodicity. There is, however, some strong hint that 
the object indeed behaves periodically in the TeV regime: the 
lower panel in Fig. Q] shows the energy flux of 1 TeV photons 
of PSRB1259-63 as a function of the separation distance r be- 
tween the pulsar and its companion star. The bidaily 2004 en- 
ergy fl uxes for 1 TeV photons are taken from lAharonian et al.l 
(2005) (see their Fig. 8) whereas the 2007 fluxes are calcu- 
lated from bolometric fluxes above 1 Te V assuming the underly - 
ing spectral power law as described in lAharonian et all ([2009 ). 
There is an apparent relation between pre and post periastron 
data obtained in the years 2004 and 2007 by H.E.S.S. and the 
orbital distance. The correlation coefficient between the 2007 
fluxes shown in time bins of months and the monthly weighted 
mean of the bidaily 2004 data as a function of r amounts to 
p = 0.89 + 0.12. Despite the limited statistics of the TeV data 
and the fact that firm conclusions can not be drawn at present, 
this behavior would indicate essentially two things: (i) the TeV 
flux is symmetric with respect to periastron and thus mainly a 
function of the binary separation, (ii) the system shows indirect 
signs of periodicity. It should be noted that in frameworks of 
a leptonic scenario this kind of behavior is rather unexpected. 
Indeed, since any two different orbital positions, which are sym- 
metric with respect to periastron, are characterized by different 
IC scattering angles (according to the currently accepted orbital 
elements), one should expect different fluxes i n the direction of 
the observer due to a nisotropic IC scattering dKirk et al.l fl999: 
lAharonian & Atovanlll98ll ). The significance of this effect de- 
pends on many different factors: the actual location of the pro- 
duction region, the y-ray energy band, the slope of the electron 
spectrum and the temperature of the optical star. Another source 
of asymmetry in PSRB 1259 -63 should be Doppler modulation 
of the post shock flow of relativistic particles. The orbital param- 
eters of the system in connection with expected large Lorentz 
factors of the bulk flow imply a strong phase dependen ce of the 
related non-thermal emission dKhangulvan et al.ll2008l) . Taking 
into account the sys tem orientation with respect to the observer 
(Johnston et a 1. 1992), this naturally should lead to an asymme- 
try in the observed VHE fluxes. In Figf2]the ratios of the VHE 
fluxes (integrated between 380 GeV and 10 TeV) for the loca- 
tions of the pulsar symmetric with respect to periastron pas- 
sage (i.e. characterized by the same absolute value of the true 
anomaly) are shown. This figure illustrates the dependence of the 
flux ratio not only on the pulsar orbital location, but also on the 
star temperature. For the calculations the electron distribution 
was assumed to be a power law with index 2. For the currently 
accepted star temperature T* = 2.3 x 10 4 K, the fluxes should dif- 
fer by a factor of 1 .5 at the locations with true anomaly 6 ~ 2rad 
(or 6 ~ 115°, r ~ 3 x 10 13 cm). Apparently, the H.E.S.S. obser- 
vations from these epochs show almost equal fluxes (see FigJTJ 
within the limited statistics of the data. In case this feature turns 
out to be significant there are several possible explanations in the 
IC scenario for such an observational result. 
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Fig. 1. Energy fluxes of 1 TeV photons from PSR B 1259-63 as a 
function of the true anomaly (top panel) and of the binary separa- 
tion distance (bottom panel). The measurements were carried out 
by H.E.S.S. during the 2004 (black full symbols - fluxes using 
a bin width of 2 days) and 2007 (red empty symbols - monthly 
fluxes) periastron passage, respectively. Squares denote pre- and 
circles post periastron measurements, respectively. The vertical 
grey boxes indicate the location of the stellar disc inferred from 
the best fit geometrical parameters from radio data. The vertical 
dashed and dotted lines denote the positions of the stellar equa- 
torial plane along t he pre- and post periastron orbit, respectively 
(Bogomaz ovl2005l) . 



(i) The ratio of the fluxes strongly depends on the temper- 
ature of the star (see Figf2|. Indeed, if the IC scattering occurs 
deeper in the Klein-Nishina regime, then the cross section shows 
a weaker angular dependence. Thus, if the star temperature were 
larger by a factor of 1.5, then the ratio of the fluxes would be 
consistent with the IC scenario, considering the relatively large 
error bars. 

(ii) If the VHE emission production site has a size not much 
smaller than the distance to the optical star, then the IC scattering 
from each location of the pulsar will be characterized by differ- 
ent scattering angles, thus the orbital phase dependence of the 
VHE flux will be smoothed. 

(iii) The effects of y-y absorption, anisotropic IC scattering 
and Doppler boosting may compensate each other. Thus the ob- 
served flux accidentally appears at the same level. It has to be 
noted, however, that as for photon absorption this possibility re- 
quires a higher luminosity of the optical s tar, otherwise the y- 
y attenuation is negligible dKirk et al.ll 19991) . 
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Fig. 2. Ratio of VHE fluxes (integrated between 380GeV- 
lOTeV) for any two different points of the orbit characterized by 
the same distance to the optical star as a function of star temper- 
ature and orbital phase. The distribution of scattering electrons 
was assumed to have a power law index equal to 2. The dashed- 
dotted line corresponds to the contour line for a flux ratio of 0.8. 
Given the relatively large statistical and systematic uncertainties 
(close to 20%) of the VHE data, in case of a star temperature 
close to (or higher than) ~ 4 x 10 4 K, the TeV light curve might 
look symmetric with respect to the periastron passage. However, 
if the star temperature is lower (e.g. ~ 2 x 10 4 K), the differ- 
ence in the fluxes for orbital locations with a true anomaly of 
6 = ±115° (or G = ±2 rad) should be measurable with H.E.S.S. 



While the second possibility requires a very detailed mod- 
elling of the system, which would cover self-consistently the 
hydrodynamics of the interaction, particle acceleration and non- 
thermal emission production, the points (i) and (iii) (except for 
Doppler boosting) may be checked through a detailed study of 
the properties of the optical star in the system. 
X-ray observations can provide important information for the 
study of the leptonic scenario. The same electrons, which emit 
VHE emission through the IC scattering, will produce X-ray 
synchrotron emission. Moreover, this channel has certain ad- 
vantages as compared to VHE emission: (i) X-rays are emit- 
ted isotropically, independent of the orbital location (unless the 
magnetic field has an ordered structure); (ii) X-rays should not 
suffer absorption (given the large size of the system). However, 
these advantages are compensated to some extent by an addi- 
tional a priori unknown parameter, the magnetic field in the 
production region. Considering the magnetic field being only 
dependent on the separation distance b etween the two objects 
dRirk et al.ll999l:lKhangulvan et al.l 2007). one should expect the 
X-ray fluxes to be symmetric with respect to periastron passage, 
as is the case for TeV y-rays. In X-rays, however, such a behav- 
ior can not b e seen. In Fig. [3] X-ray data in the 1 — 10 k eV band 
from ASCA (iKaspi et al.f 1 1 994t iHiravama et al]|j996l) Beppo- 
Sax, XMM Newton dChernvakova et alj 120061 120091). Swift 
Chan dra and Suzaku (Uchivama et al. 2009; IChernvakova et al.l 
2009) from the last periastron passages since 1994 are shown 
as a function of the binary separation r. Here the picture is dif- 
ferent from the symmetric lightcurve seen in the VHE regime. 
For r > 2 x 10 13 cm pre- and post periastron phases are sep- 
arated by almost one order of magnitude in flux since there is 
a steep drop in the pre periastron part of the lightcurve (black 
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Fig. 3. X-ray fluxes in the 1-10 keV band from pre periastron 
(black symbols ) and post periastron (red symbols ) phases mea- 
sured in the years 1994, 1997, 2001, 2004 and 2007 as a function 
of the binary separation distance r. Also shown are three differ- 
ent model X-ray light curves corresponding to the adopted cool- 
ing profiles shown in Fig. [4] The grey band indicates the exten- 
sion of the stellar disc. The vertical dashed and dotted lines cor- 
respond to the location of the stellar equat orial plane along t he 
pre- and post periastron orbit, respectively (Bog omazovl2 005). 



symbols ). Interestingly, this point of broken symmetry in the X- 
ray lightcurve is roughly coincident with an increase in the TeV 
flux by a factor of ~ 3 (compare Fig. [I). The clear discrepancies 
between the X-ray and the y-ray lightcurves presented here may 
indicate that the non-thermal particles responsible for the corre- 
sponding emission could be different or that the magnetic field 
has a rather complicated dependence on the orbital phase. 



3. Description of the model 

In this section, the model that was used to calculate the TeV en- 
ergy flux of PSRB 1259-63 taking into account adiabatic cool- 
ing of VHE el ectrons is briefly explain ed. For a more detailed 
description see Khangu rvan et al.l (120071) . 

The electron steady state distribution function (EDF) n e (t, y) 
of shock accelerated electrons as a function of time t (i.e. or- 
bital position) and energy y can be described by t he following 
expression (see e.g. Ginzburg & Syrovatskii ( 1964)): 



n e (t, y) 



1 



\y(y, 01 



oo 

J Qiy'W 



(i) 



Here, the electron energy loss rate y due to e.g. synchrotron ra- 
diation (syn), IC cooling (ic) and adiabatic losses (ad) can be 
written as 



y = Tsyn + tic + tad • 



(2) 



While th e calculation of y svn oc B 2 y 2 and y lc i s straight forward 
(see e.g. iBosch-Ramon & KMng ulvanl d2009l) ) it is highly non 
trivial to compute tad analytically since it depends on the regime 
of the hydrodynamical interaction between the PW and the stel- 
lar outflow. Therefore, the adiabatic cooling profile is inferred 
by comparing an interpolation of the VHE data with the ex- 
pecte d lightcurve for c onstan t electron injection, i.e. y a d = 
(see lKhangulvan et all d2007h their figure 2). 



The other notations in Eq. ([TJ are the following. The source 
term Q(y) denotes the incident injection spectrum of shocked 
electrons from the pulsar wind. Given the lack of understanding 
of the acceleration processes in binary systems, it was assumed 
to be a simple power law with an energy cutoff at E efma : 



Q{t,y)=Ay- a e 



Ky-y mm) 



(3) 



where 9 is a Heaviside function. In general, this kind of accel- 
eration spec trum with a ~ 2 is consi stent with observations of 
pulsars (e.g. lKennel & Coronitil[l9 84). The spectrum is defined 
by the following parameters: A is the normalization coefficient, 
determined by the fraction of PW particles transferred to VHE 
injection particles; a is the slope of the accelerated spectrum 
and is defined by the acceleration mechanism, and y m j n is the 
minimum energy in the accelerated spectrum. The cutoff energy 
£e,max should be determined by the balance between acceleration 
and losses. A simple estimate for the corresponding acceleration 
timescale, which is a function of the magnetic field B(r), which 
itself depends on the binary separation r as B(r) oc 1 jr , is 



Rl 



(4) 



where = mcyjeB(f) is the Larmor radius for leptons and 
r\ > 1 characterizes the acceleration rate. Comparison of Eq. (f2]) 
and Eq. (0]i provides the value of E e>max , The EDF obtained with 
the inferred adiabatic loss profile was used to calculate the IC 
fluxes. Since adiabatic losses are assumed to be dominant, the 
EDF is nearly independent of the strength of the magnetic field. 
Attenuation of the TeV emission due to y -y absorption ap pears 
to be non-significant for PSRB 1259-63 dKirk et alj|1999l) and 
thus has not been considered in the calculation for the VHE 
fluxes presented here. 

4. Results 

In this study, three different time profiles for the non radiative 
cooling coefficients as shown in Fig. [4] were adopted in order to 
account for the y-ray flux variability in PSRB1259-63 as mea- 
sured by H.E.S.S. For the calculations, the following orbital pa- 
rameters were used: eccentricity e = 0.87, periastron separation 
ro = 9.6 x 10 12 cm, apastron s eparation r a = 1 .4 x 10 14 cm and in- 
clination of the orbit i 35° dJohnston et alJ[l992l) . Using these 
non-radiative loss profiles for y a d = y a d(0 in Eq. d2J, the EDF 
n e (t,y) is obtained through Eq. ([]]). Thus, calculations of the IC 
energy flux using the newly computed EDFs that additionally 
account for adiabatic losses give the predictions shown in Fig. [5] 
For this calculation, it was assumed that 10 % of the pulsar spin- 
down luminosity was transferred to non-thermal particles. The 
magnetic field at periastron was chosen to be B = 0.1 G. An 
energy range of 1 GeV to 10 TeV for the injected electron spec- 
trum was used. Other p a ramet ers were assumed to be the same 
as in Khangul van et al.l d2007h . In comparison, the latter study 
employed an adiabatic loss profile inferred from H.E.S.S. 2004 
data only, not accounting for the drop in the VHE photon flux 
at per iod r « 17 d corresponding to 9 a +90°. Uchivama et al.l 
(2009), on the other hand, assumed the time scale of adiabatic 
losses around r ~ 30 d to be ~ 10 times shorter than indicated 
in this study in order to match the broadband spectral energy 
distribution of the system while otherwise adopting similar pa- 
rameters and radiation processes. 

Comparing Fig.[T](bottom panel) and Fig.|5]it seems that there is 
a peculiar "dip" in the TeV emission around r ~ 2 x 10 13 cm that 
is best reproduced when assuming a non-monotonic loss rate as 
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Fig. 4. Cooling coefficients for PSRB 1259-63 as inferred from 
the combined H.E.S.S. datasets shown in Fig.Q] The grey bands 
indicate the extension of the stellar disc. The vertical dashed and 
dotted lines correspond t o the location of the stellar equatorial 
plane (Bogomazov 2005). 
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Fig. 5. Predicted energy fluxes of 1 TeV photons from 
PSRB 1259-63 corresponding to the three different non radia- 
tive loss profiles as shown in Fig.|4]compared to the 2004 (black 
dots) and 2007 (red empty squares) H.E.S.S. data. The vertical 
grey boxes indicate the location of the stellar disc. The vertical 
dashed an d dotted lines deno te the position of the stellar equato- 
rial plane ( Bo gomaz ovl2005h . 

depicted by the black solid curve in Fig. |4] Using a mono tonic 
and thus less complicated cooling function such as the green 
dash dotted profile gives also a good overall agreement with the 
H.E.S.S. data but fails to describe the reduced photon fluxes seen 
at 9 « ±90°. 

The predictions of X-ray fluxes (1-10 keV) for the adopted non- 
radiative loss profiles and model parameters are shown in Fig. [3] 
The pre periastron data (black symbols) deviate considerably 
from the model, which adopts a symmetric behavior of the B- 
field with respect to periastron. Obviously this assumption alone 
does not explain the X-ray data. 

5. Discussion & Conclusions 

PSRB 1259-63 has been observed for several orbits over the 
past decades exhibiting a globally consistent picture for the 



phase dependent electromagn etic radiation at radio wavelengths, 
X-Rays and VHE (see e.g. lUchivama et all (l2009t) their Fig. 2). 
This suggests a physical link between the emission mechanisms 
in the different energy regimes such as the increasing radiation 
efficiency with decreasing binary separation distance r resulting 
in enhanced PW and stellar wind interactions. However, a more 
detailed comparison between the VHE and X-Ray lightcurves 
as a function of r reveals substantial qualitative differences in 
terms of symmetry with respect to periastron and evolution in 
flux level. While in the X-Ray lightcurve, Fig. [3] pre- and post 
periastron data are separated by almost one order of magnitude 
in flux level for r > 2 x 10 13 cm such an asymmetry can not 
be seen in the TeV band as shown in Fig. [T] Even if the TeV 
data suffers from low statistics, a similarly prominent feature in 
this energy regime should be notable. This is also reflected in 
the unexpected high TeV flux at r ~ 3.6 x 10 13 cm (6 ~ -125 
deg) when a ssociating the TeV emissi on with the disc location 
proposed in Chernvakova et al. (2006). The latter was inferred 
mainly from post periastron data (see lAharonian et al ] (I2009h 
their Fig. 7) demanding for an asymmetric lightcurve due to a 
tilted stellar disc with respect to the pulsar orbital semi minor 
axis. However, it has to be noted that at present the discrepancy 
of the X-ray and TeV data, in the above discussed context, 
relies mainly on this single pre periastron VHE point. Another 
deviation in the two energy bands relates to the increase in 
pre periastron X-ray fluxes seen for decreasing r in the range 
2.4 x 10 13 cm > r > 1.5 X 10 13 cm. This evolution appears to 
be reversed in the VHE lightcurve where there is a significant 
drop in the TeV photon flux between r ~ 2 x 10 13 cm and 
r ~ 1.5 x 10 13 cm. 

Taking into account the orbital orientation of 
PSRB 1259-63 with respect to the observer, Doppler boosting 
of post shock particle flows should affect the TeV lightcurve and 
yield strong asymmetries which are not seen in Fig.Q] Of course 
Doppler boosting could in principle always be accidentally 
compensated by other factors such as anisotropic IC. In general, 
however (and especially in case of an IC scenario), the VHE 
flux from a binary system is not expected to be symmetric with 
respect to the periastron passage. Indeed, the corresponding 
symmetric orbital phases use to be characterized by different 
y-y opacities and scattering angles with respect to the line of 
sight. The flux symmetry with respect to periastron passage, 
indicated in the case of PSRB 1259-63, suggests either a 
rather fine compensation of the above mentioned effects or 
(i) a negligibly small absorption of the y-rays and (ii) a y-ray 
production mechanism, which has a weak dependence on the 
orbital phase. The latter could be the case for e.g. a large 
production region in PSR B1259-63 or IC scattering to proceed 
in the deep Klein-Nishina regime with the cross section a- K only 
changing marginally with 9. 

To explain the variability of the observed flux as a func- 
tion of the separation distance one possibility is to intro- 
duce non-radiative (adiabatic) losses, which would dominate 
in PSRB1259-63 over the whole orbital period (see Fig. 4 in 
KhangulyanetaD (120071) ). Moreover, the observed lightcurve 
with a number of humps and dips suggests a rather compli- 
cated dependence of the non-radiative losses on the separation 
distance. Given the complexity of the self-consistent calculation 
of the adiabatic losses, TeV data were used to infer a possible 
profile of the losses. In particular, one may expect the adiabatic 
loss rate to have a peak close to periastron together with two 
smaller peaks located at orbital positions characterized by a true 
anomaly of 6 « +75°. Those smaller peaks may be linked to the 
impact of the stellar disc, as indicated in Fig. |4] when the pulsar 
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exits the equatorial wind and thus the loss rate due to interfer- 
ence with the outflow decreases. 

The predicted TeV lightcurves including different possible cool- 
ing profiles, Fig.Hl give an overall qualitative agreement with the 
data (see Fig. [5j>. Profiles with a simple evolution of the cooling 
rate such as the dashed curves in Figs. [4] and [5] do not account for 
the TeV data close to periastron (9 * +70°). Naturally, the best 
agreement is achieved by a cooling function featuring two ad- 
ditional peaks that account for the potential impact of the stellar 
disc (black solid curve in Figs.|4]and[5]i. The predicted lightcurve 
shows the moderate impact of anisotropic IC scattering while 
still being qualitatively compatible with the observations. 
Regarding X-ray emission, the predicted lightcurves correspond- 
ing to the presented cooling coefficients only show a weak de- 
gree of agreement to the data (Fig. [3}. The most prominent dis- 
agreement stems from pre periastron data for binary separation 
distances r > 2 x 10 13 cm. This basically indicates a more com- 
plicated r-dependence of the magnetic field in the production re- 
gion than the symmetric B(r) oc 1/r assumed in this study. There 
is qualitative agreement for the post periastron lightcurve (red 
symbols) as far as flux level and global evolution is concerned. 
Introducing a second peak close to periastron in the black solid 
model curve seems also not justified in X-rays even if improv- 
ing the prediction compared to the red dashed simple curve. In 
conclusion, none of the suggested model curves accounts quan- 
titatively for the observational data in this energy regime based 
on a simple phase dependence of the B-field. This requires addi- 
tional assumptions on the orbital phase dependence of the mag- 
netic field. 

Although any detailed discussion of the ratio of TeV and X-ray 
fluxes requires rather accurate calculations and goes beyond the 
scope of this paper, a qualitative explanation for a sharp increase 
of the X-ray flux before periastron passage may be suggested 
here. In the framework of dominant non-radiative losses, a de- 
crease of the IC flux by a factor of 3 should be caused by an 
equivalent enhancement of the non-radiative losses. This may be 
achieved with an increase of the ram pressure in the stellar out- 
flow, so that the PW termination shock moves closer to the pulsar 
(roughly by a factor of 3). Such a scenario is naturally provided 
by the pulsar crossing the dense stellar disc. This would not sig- 
nificantly change the density of the target photons (assuming the 
shock is located close to the pulsar), but should lead to a sig- 
nificant increase of the magnetic field in the production region 
(again by a factor of 3). Thus, it is indeed natural to expect an 
increase of the synchrotron flux by one order of magnitude. On 
the other hand, it has to be noted that the overall orbital behavior 
of the X-ray flux cannot be explained by means of these sim- 
ple arguments. In this regard it is noteworthy that the increase in 
pre perisatron X-ray emission coincides with the position of the 
stellar disc (see black points and vertical dashed line in Fig. [3]). 
The steepest increase of the emission is roughly aligned with the 
equatorial stellar plane where the disc density is presumed to be 
highest. Again, it is difficult to account for the behavior of the 
post periastron data (red points) in this model picture. 

Although the H.E.S.S. data are not yet significant enough 
to be conclusive regarding the existence of two dips in the TeV 
lightcurve, the symmetry of the presented feature in the VHE 
regime with respect to periastron and the correlation with a 
significant rise of the X-ray emission at the same separation 
distance are certainly notable. It could be well explained by a 
non monotonic adiabatic cooling profile, tracing a change in the 
shock region's size and magnetic field conditions induced e.g. by 
stellar matter outflows of increased density such as the stellar 
disc. 



Future observations in the VHE regime with instruments 
such as CTA or H.E.S.S. II should shed light on this interesting 
question. 
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